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Abstract. Trying to obtain a more complete picture of 
star forming regions in the Taurus molecular cloud, four 
newly discovered dense molecular cores (L1521D, L1521F, 
L1524, L1507A) are identified and mapped through the 
ammonia (J,K) = (1,1) and (2,2) rotational inversion lines. 
These cores have sizes from 0.06 to 0.09 pc, hydrogen den¬ 
sities from 0.6 10 4 to 19.9 10 4 cm” 3 and kinetic tempera¬ 
tures between 7.9 and 9.9 K. The masses range from 0.2 
and 1.0 Mq, placing the cores with the lowest M values at 
the lower edge of the mass distribution for ammonia cores 
in the TMC. Turbulent, thermal and gravitational ener¬ 
gies have been estimated. A comparison between these en¬ 
ergy terms and considerations related to thermal and tur¬ 
bulent line broadening indicate that these cores are close 
to gravitational equilibrium. Moreover, we report detec¬ 
tions of the J=9-8 transition of HC5N towards the am¬ 
monia peak positions of these four molecular cores. The 
HC5N column densities are between 1.6 10 12 and 9.2 10 12 
cm” 2 , in agreement with values derived for other molecu¬ 
lar cores located in the TMC. 

Key words: ISM: clouds - ISM: molecules - Radio lines: 
ISM - ISM: individual objects: TMC 


1. Introduction 


from the IRAS Point Source Catalogue (PSC, IRAS 1985), 
molecular outflows (e.g. Fukui 1989| ; Fukui et al. 1993) and 
water maser sources (Wouterloot et al. 1993| ), indicating 
ongoing star formation. The TMC contains several dark 
molecular clouds of siz e ~ 1 pc and densities of 10 3 cm” 3 
(e.g. Walmsley et al. 198 ()t W outerloot & Habing 1985 


Ungerechts & Thaddeus 1987 ; Mizuno et al. 1995| ). These 
in turn contain dense molecular cores, high density (> 10 4 
cm” 3 ) subs tructures which have b een cl assified by Myers 
et al. ( |1983 ) and Myers & Benson ( 1983 ). For the sky dis- 
tribut ion of the molecular cores in the TMC, see Kenyon 
et al. ( 199C ). Sub sequen t mult imole cular observation s (e.g . 
Benson & Myers 1984: Fuller 198S; Fuller & Myers 1993 ) 
have confirmed that dense cores are smaller than 0.1 pc 
with masses ranging from 0.5 to 10 solar masses and ki¬ 
netic temperatures of the order of 10 K. Beichman et al. 
( [1984 ; 1986 ) have shown that about 50% of the molecular 
cores are associated with IRAS sources, indicating that 
these objects may contain, in their interiors, young stellar 
objects (YSOs). Therefore, it is now clear that there exists 
a close association between low-mass YSOs and molecu¬ 
lar clumps. Since the TMC is not strongly affected by 
the nearby Cas-Tau OB association (Blaauw 1991), it is 
clear that it provides unique opportunities to study spon¬ 
taneous low-mass star formation and, in particular, to 
investigate the nature of the earliest stages of low-mass 
star formation, strictly connected with dense molecular 
clumps. 


The Taurus molecular cloud (TMC) is one of the best 
targets to study low-mass star formation, since it is lo¬ 
cated at a distance of only 140 pc. The region is asso¬ 
ciated with T Tauri stars, identified by infrared and op¬ 
tical observations (e.g. Strom et al. 198£; Kenyon et al. 
199C; Weintraub 1990|), low-luminosity objects extracted 
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Following the definition given by Myers & Benson 


(1983; i.e. sources with T\ > 0.4 K at the Haystack and 
Green Bank antennas), 16 dense molecular cores have so 


far been identified in the TMC (Myers & Benson 1983 
Benson & Myers 1989|). In order to obtain a more com¬ 


plete sample of TMC cores, to gain information about 
the nature of the molecular gas before or at the onset of 
YSO formation and to clarify the evolutionary stages con- 
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nected with the collapse process, we report observations 
of the NH 3 (J, K) = (1,1) and (2,2) inversion lines and the 
cyanodiacetylene (HC 5 N) J=9-8 rotational line towards 4 
newly discovered dense cores. 


2. Selection criteria 


The molecular cores chosen for this project have been se¬ 
lected from an extensive search for ammonia emission to¬ 
wards 26 high visual extinction regions in the TMC iden¬ 
tified from the Palomar Observatory Sky Survey (POSS) 
maps and from Myers et al. ( 197S| ). The observed posi¬ 
tions are located within 04 /l 15 m < a < 05^ 00 m and 
+16° < 5 < +30°: 9 regions in the long filament which 
includes TMC-2 (Kutner’s Cloud), 7 regions in the paral¬ 
lel filament to the north, containing B217 (see Fig. 1 of 
Walmsley et al. 1980D , and 10 isolated optically obscured 
regions where no extended clouds were seen on the POSS 
plates. For the survey, the Haystack (USA) 37-m and the 
Green Bank (USA) NRAO 43-rn antennas have been used 
with FWHP (full width half power) beam sizes of 84" and 
beam efficiencies of about 0.25 during several runs in 1983. 

Table 1 summarizes the regions where NH 3 line emis¬ 
sion was searched for. For the non-detections the r.m.s. 
T'a value is reported. Ammonia emission has been de¬ 
tected towards 9 sources. Out of these 9, 5 can be classified 
as dense molecular cores: L1521D, L1521F, L1524, L1535 
and L1507A. Since L1535 has already been investigated 
by Ungerechts et al. ( 1982| ), we have focused our atten¬ 
tion on the other four cores, which have been mapped in 
NH 3 , using the 100-m MPIfR radio-telescope at Effelsberg 
(Germany). 


Table 1. Regions investigated through NH 3 observations 


a 

(1950) 

S 

(1950) 

n 

(mK) 

Antenna ^ 

Dense 

core 

04 

17 

56 

+26 

55 

23 


449 

GB 

L 

1521 

D 

04 

18 

55 

+29 

21 

00 

< 

95 

HA 




04 

19 

55 

+19 

19 

49 

< 

102 

GB 




04 

20 

30 

+26 

32 

47 

< 

101 

GB 




04 

21 

48 

+27 

10 

58 

< 

95 

GB 




04 

23 

36 

+24 

32 

40 

< 

150 

HA 




04 

25 

00 

+18 

35 

12 

< 

103 

GB 




04 

25 

35 

+26 

45 

00 


977 

GB 

L 

1521 

F 

04 

26 

20 

+24 

28 

06 


603 

GB 

L 

1524 


04 

26 

25 

+26 

52 

47 

< 

102 

GB 




04 

27 

58 

+24 

21 

06 

< 

113 

HA 




04 

28 

02 

+24 

05 

20 

< 

108 

HA 




04 

28 

40 

+27 

02 

13 

< 

85 

GB 




04 

29 

38 

+24 

42 

13 

< 

80 

HA 




04 

29 

40 

+24 

06 

00 

< 

151 

HA 




04 

29 

57 

+24 

10 

53 

< 

179 

HA 




04 

31 

08 

+24 

02 

10 


216 

GB 




04 

32 

37 

+24 

01 

42 


788 

HA 

L 

1535 


04 

37 

28 

+29 

50 

00 

< 

74 

HA 




04 

37 

39 

+29 

48 

37 


270 

GB 




04 

39 

29 

+29 

38 

07 


472 

GB 

L 

1507 

A 

04 

40 

20 

+29 

34 

37 


276 

GB 




04 

42 

18 

+16 

31 

58 

< 

85 

GB 




04 

42 

22 

+17 

01 

42 

< 

86 

GB 




04 

42 

38 

+16 

51 

35 

< 

85 

GB 




04 

57 

16 

+26 

10 

23 


236 

GB 





^ Observed with the 43-m Green Bank antenna (GB) or with 
the 37-m Haystack antenna (HA) in 1983. The FWHP beam 
sizes are ~ 84 arcsec, while the beam efficiencies are ~ 0.25. 


3. Observations 

The observations were made during several runs in 1985, 
1986 and 1996 with the MPIfR antenna employing a 
maser receiver during the 1985-1986 and a HEMT re¬ 
ceiver during the 1996 periods. At the rest frequency of 
the (J, A") = (1,1) and (2,2) lines of NH 3 (23694.496 MHz 
and 23722.631 MHz, respectively) and of the J=9-8 line 
of HC 5 N (23963.888 MHz) the FWHP beam size is 40". 
The zenith system temperature was around 60 K in good 
weather conditions on an atmosphere corrected antenna 
temperature (T£) scale. The observations were made using 
position switching, with an integration time of typically 
3 minutes on source. Spectral information was obtained 
with a 1024-channel three-level autocorrelation spectrom¬ 
eter. We observed the NH 3 (1,1) and (2,2) lines simulta¬ 
neously by splitting the autocorrelator into two spectrom¬ 
eters of 512 channels with bandwidths of 6.25 MHz and 
channel spacings of 0.154 km s -1 . The velocity resolution 
was 0.185 km s^ 1 , i.e. a factor 1.2 larger than the channel 
spacing. We observed the HC 5 N line with the autocorre¬ 
lator split into two parts, each with a bandwidth of 3.12 
MHz and a velocity resolution of 0.092 km s -1 . The spec¬ 


tra from the two split parts were then co-added in order 
to reduce the noise. The intensity scale of the spectra was 
calibrated on the continuum sources NGC7027, 3C123 and 
3C2 86, ac cording to the flux densities given by Baars et 
al. ( 1977 ). The pointing accuracy was better than 10", 
corresponding to a point source flux density error of < 
15%. The calibration accounted for the dependence of the 
gain on the elevation. Total flux density uncertainties are 
of the order of 20 %. 


4. Results 

Ammonia maps have been obtained of the four dense cores 
L1521D, L1521F, L1524 and L1507A in the (J, I <) = (1,1) 
line (Figs. 1 to 4). In order to obtain the observed source 
size, we have determined the half maximum contours of 
the ammonia main beam brightness temperature (Tmb! 
K) and have calculated the apparent geometric mean di¬ 
ameter (FWHP) D = y/ab : where a and b are the major 
and minor-axis, respectively. The ammonia core sizes are 
0.09 pc (L1521D), 0.07 pc (L1521F), 0.08 pc (L1524) and 
0.06 pc (L1507A). Deconvolution by the 0.027 pc beam 



















Codella et al.: Four new dense molecular cores in the TMC 


3 


Table 2. Observed NH 3 (1,1) line parameters. The (1950) coordinates of the peak positions of the four cores are: (a = 04 fe 
17™ 59 s ; <5 = +26° 55' 43") and (a = 04 h 18 m 02 s ; 8 = +26° 56' 03") for L1521D, (a = 04 h 25™ 36 s ; <5 = +26° 45' 00") for 
L1521F, (a = 04 h 26™ 21 s ; 8 = +24° 28' 26") for L1524 and (a = 04 h 39™ 29 s ; 8 = +29° 38' 17") for L1507A 


Name 

Aa 

(") 

A<5 

(") 

Tmb 

(K) 

ULSR 
(km s” 1 ) 

FWHM 
(km s” 1 ) 

Ttot 

Avint 
(km s” 1 ) 

L1521D 

+40, 

+20 

2.3 (0.1) 

6.94 (0.15) 

0.8 (0.2) 

4.5 ( 1 . 1 ) 

0.33 (0.02) 


+80, 

+40 

2.3 (0.1) 

6.97 (0.15) 

0.6 (0.2) 

4.6 (1.3) 

0.30 (0.02) 

L1521F 

+20, 

0 

6.1 (0.3) 

6.62 (0.15) 

0.7 (0.2) 

12.0 (0.8) 

0.25 (0.01) 

L1524 

+20, 

+20 

4.7 (0.2) 

6.40 (0.15) 

0.8 (0.2) 

14.0 (0.9) 

0.26 (0.01) 

L1507A 

0 , 

+ 10 

3.4 (0.2) 

6.16 (0.15) 

0.7 (0.2) 

5.8 (0.9) 

0.27 (0.01) 


L S2TD 



CD 

OT 

Fig. 1. Contour NH 3 (1,1) Tmb map of L1521D. The (0,0) <j 
position corresponds to: a = 04 h 17™ 56 s ; <5 = +26° 55' 23". 

The contour levels are from 1.4 K to 2.4 K by steps of 0.2 K 
(~ 1 r.m.s. noise). The empty circle shows the Effelsberg beam 
(HPBW), the small crosses mark measured positions. 


does not significantly alter these values. For those spectra 
where the main and satellite hyperfine components have 
been observed with a sufficient signal-to-noise ratio (> 3) 
we have obtained the total optical depth (T tot ) and the 
intrinsic linewidth (FWHM) of the hyperfine components 
(Avint; km s^ 1 ) using a program which fits these param¬ 
eters to the observed spectra, assuming gaussian profiles 
in an individual component and relative opacities being 
consistent with the Local Thermodynamic Equilibrium 
(LTE) approximation. For the positions where the satel¬ 
lite components have intensities below the noise level, but 
the main component of the NH 3 (1,1) line has been de¬ 
tected, we have fitted a single gaussian profile. Moreover, 
we have observed the NH 3 (2,2) transition towards the 
peak positions of the four molecular cores. 


Table 3. Observed NH 3 (2,2) line parameters 


Name 

Aa 

(") 

A<5 

(") 

Tmb 

(K) 

VLSR 

(km s -1 ) 

FWHM 
(km s” 1 ) 

L1521D 

+40, 

+20 

0.5 (0.1) 

6.81 (0.15) 

0.6 (0.2) 

L1521F 

+20, 

0 

1.3 (0.2) 

6.65 (0.15) 

0.3 (0.2) 

L1524 

+20, 

+20 

0.9 (0.2) 

6.44 (0.15) 

0.3 (0.2) 

L1507A 

0 , 

+10 

1.1 (0.3) 

6.17 (0.15) 

0.1 (0.2) 


LI 521F 



100 50 0 -50 - 100 

R.A. Offset (arcsec) 


Fig. 2. Contour NH 3 (1,1) Tmb map of L1521F. The (0,0) 
position corresponds to: a = 04 h 25™ 35 s ; 8 = +26° 45' 00". 
The contour levels are from 2.6 K to 6.4 K by steps of 0.4 I< 
(~ 2 r.m.s. noise). The empty circle shows the Effelsberg beam 
(HPBW). 


Table 2 displays the observed NH 3 (1,1) line param¬ 
eters regarding the peak positions: the source name, the 
right ascension and declination offset (") of the peak po¬ 
sition with respect to the coordinates given in Table 1, 
the main beam brightness temperature, the LSR velocity 
(ulsr; km s- 1 ), the apparent line width (FWHM; km s” 1 ) 
of the main component, the total optical depth and the 
intrinsic linewidth of the hyperfine components including 
instrumental broadening. The errors (r.m.s.) of the de¬ 
rived parameters are given in parentheses. Figure 5 shows 
the NH 3 (1,1) and (2,2) spectra of the peak position of 
each mapped source. Table 3 lists the observed NH 3 (2,2) 
gaussian line parameters towards the NH 3 (1,1) peak po¬ 
sitions: the main beam brightness temperature, the LSR 
velocity and the FWHM line width including instrumental 
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broadening. For L1521D, we report the observed parame¬ 
ters towards the (+40, +20) peak position, since the inte¬ 
grated ammonia emission associated with this position is 
slightly larger that that of the (+80, +40) peak position. 



100 50 0 -50 

R.A. Offset (arcsec) 


Fig. 3. Contour NH3 (1,1) Tmb map of L1524. The (0,0) posi¬ 
tion corresponds to: a = 04 h 26 m 20 s ; 5 = +24° 28' 06". The 
contour levels are from 2.1 K to 4.5 K by steps of 0.4 K (~ 
2 r.m.s. noise). The empty circle shows the Effelsberg beam 
(HPBW). 


Table 4 gives the HC 5 N (J=9-8) line parameters to¬ 
wards the ammonia peak positions, obtained from gaus- 
sian fits: the parameters are main beam brightness tem¬ 
perature, LSR velocity and FWHM line width including 
instrumental broadening. Figure 6 shows the HC 5 N (J=9- 
8 ) spectra. 


5. Discussion 


5.1. Derived parameters from NH 3 observations 

The rotational temperature (T ro t) is defined by applying 
the Boltzmann distribution to the populations of different 
(J,K) rotational levels and has been calculated from the 
( 1 , 1 ) and ( 2 , 2 ) data follo wing the proc edure outlined by 
Hiittemeister et al. ( 1993 ) and Lemme ( 1995 ). Given the 
rotational temperatures, it is possible to estimate the ki¬ 
netic temperature (Tkin)- Since the derived T rot values (see 
Table 5) are all around 10 K, the differences between ki¬ 
netic and rotational temperatures, as shown by Walmsley 
& Ungerechts (1983) and Danby et al. (1988), are within 


L1507A 



R.A. Offset (arcsec) 


Fig. 4. Contour NH3 (1,1) Tmb map of L1507A. The (0,0) 
position corresponds to: a = 04 h 39 m 29 s ; 8 = +29° 38' 07". 
The contour levels are from 2.0 K to 3.4 K by steps of 0.2 K 
(~ 1 r.m.s. noise). The empty circle shows the Effelsberg beam 
(HPBW). 


the errors of our T rot measurements. We thus can assume 

Tjcin = Trot- 

The excitation temperature (T ex ), derived assuming a 
beam-filling factor of one, is defined in an analogous way to 
T r ot, but it is related to the population across an inversion 
doublet (J,K). Using T ex and T ro t, it is possible to derive 
the total ammonia column density (Wot; cm -2 ) following 
the equations given by Ungerechts et al. (1986). For this, 
we have neglected non-metastable (J > K) inversion dou¬ 
blets and have assumed a common rotational temperature 
for all the metastable (J = K) states. These are reasonable 
approximations for rotational temperatures around 10 K 
(e.g. Harju et al. 1993 ). The derived temperatures and col¬ 
umn density values (see Table 5) are consistent with NH 3 
observations towards o ther molecular cores in the T MC 
(e.g. Myers & Benson 1983; Benson & Myers 1989) and 


more isolated Bok globules (Lemme et al. |l996|; Scappini 


& Codella 1996). The measured (1,1) excitation tempera¬ 
tures are only slightly smaller than the rotational temper¬ 
atures (see Table 5) suggesting that the (1,1) transition is 
close to being thermalised. T ex ~ T rot also requires that 
the beam-filling factor is near unity. This implies that, for 
the sources of the present sample, most of the ammonia 
flux is arising from an extended structure. 

By balancing collisions and stimulated emission 
against spontaneous emission, the hydrogen density (tih 2 i 
cm -3 ) can be estimated (Ho & Townes 1983). The derived 
7 Th 2 values, also accounting for photon trapping, range 
from 0.6 10 4 to 19.9 10 4 cm -3 (see Table 6 ) and confirm 
that ammonia traces regions with densities > 10 4 cm~ 3 . 
Following Ho & Townes ( 1983 ), the errors on «h 2 are 17% 
and 64% for L1521D and L1507A, respectively. It is worth 
noting that the errors associated with tih 2 become large 
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Fig. 5. Ammonia (1,1) an d (2,2) spectra of the peak position of each mapped core. Source name and angular offset map position 
(in arcseconds) are indicated. 


Table 4. Observed and derived HC 5 N (J=9-8) line parameters 


Name 

Aa 

(") 

A 5 
(") 

7mb 

(K) 

ULSR 
(km s” 1 ) 

FWHM 
(km s" 1 ) 

NhcsN 
( 10 12 cm -2 ) 

L1521D 

+40, 

+20 


6.86 (0.09) 

0.4 (0.1) 

4.5 (0.7) 

L1521F 

+ 20 , 

0 


6.22 (0.09) 

0.5 (0.1) 

9.2 (1.0) 

L1524 

+ 20 , 

+20 

1.3 (0.2) 

6.16 (0.09) 

0.4 (0.1) 

7.5 (0.9) 

L1507A 

0 , 

+10 

0.6 ( 0 . 1 ) 

6.03 (0.09) 

0.3 (0.1) 

1.6 (0.3) 
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Fig. 5. (continued) 


when T ex is close to T Iot . Thus, the hydrogen densities 
derived for L1521F and L1524, i.e. for those sources with 
T ex ~ T rot , may be lower limits only. 

It is well known that intrinsic motions in molecular 
clouds and, in particular, in molecular cores are complex. 
The spectral lines of thermally excited transitions also 
contain a non-thermal contribution to the linewidth (e.g. 
Myers et al. 1991 ), which can provide information about 
the dynamics of the molecular gas. In order to derive Av, 
the actual intrinsic linewidth of the ( 1 , 1 ) line, we have 
corrected the observed intrinsic linewidth AriiNT for the 
spectral resolution of the autocorrelator used during the 


observations, Au aut , by means of Av 2 = Awj NT — AL( ut , 
assuming the filter (and the line) being gaussian shaped. 
The actual intrinsic linewidth Av is composed of a ther¬ 
mal component (A-ut) and a non-ther mal one (AvntL 
Av 2 = Av\ -I- A«^ t (e.g. Myers et al. 1991 ). It is possi¬ 
ble to see (Tables 2 and 5) that: (i) the actual intrinsic 
linewidths Av are between 0.17 and 0.27 km s -1 (Acint 
~ 0.28 km s _1 ) and (ii) the non-thermal velocity com¬ 
ponent is small, with a ratio between Av and the ther¬ 
mal linewidth component Av/Avt — 1.3. The results 
of t he NH 3 observations, performed by Benson & My¬ 
ers (198S|) using the larger beamwidths of the Haystack 
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Fig. 6. Cyanodiacetylene (J=9-8) spectra at the ammonia peak position of each molecular core. Source name and angular 
offset (in arcseconds) relative to the (0,0) positions in the ammonia maps (see Figs. 1-4) are indicated. 


and Green Bank antennas, indicate that, on average, am¬ 
monia cores with embedded IRAS source(s) have larger 
linewidths (0.45 km s -1 , not corrected for the spectral res¬ 
olution) than cores without IRAS point source (0.27 km 
s _1 ). Moreover, Myers et al. (1991) studied the linewidths 
of 61 dense condensations associated with star forming re¬ 
gions. They showed (see their Fig. 1) that the non-thermal 
component of the NH 3 core motions increases with IRAS 
luminosity more rapidly than does the thermal component 
and that in cores which can form stars more massive than 


2 Mq, non-thermal motions dominate thermal motions 
( Av/Avt ~ 4). Thus, the linewidths derived suggest that 
our NH 3 clumps are either not associated with YSOs or 
that they are sites of low-mass star formation and are pre¬ 
sumably associated with IRAS poin t sour ces of low lumi¬ 
nosity (Lfir < 10 A©; Myers et al. 1991 ), corresponding 
to a stellar mass of less than 2 Mq . 


The LSR velocity information of the four mapped cores 
has been analysed in order to search for velocity gradients: 
ulsr does not vary greatly within L1524 and L1507A, 
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since the maximum differences between the wlsr values 
of the map positions are 0.16 km s -1 and 0.14 km s , 
respectively. These values are comparable to the channel 
spacing of 0.154 km s -1 . The map of L1521F shows a sig¬ 
nificant variation (0.46 km s -1 ), but there is no clear indi¬ 
cation for a systematic velocity gradient along a preferred 
axis. In contrast, L1521D displays a large variation of ulsr 
across the map revealing a velocity gradient roughly along 
the north-south axis with a change of about 0.43 km s -1 
over 120" (cf. Fig. 1; weak line emission from outside the 
outer 1.4K contour supports this trend). This corresponds 
to ~ 5 km s -1 pc , which is a value definitely larger than 
those reported by Goodman et al. ( |1993 ) for dense cores 
(< 4 km s -1 pc -1 ). It is also worth noting that the di¬ 
rection of the gradient is not related to the major-axis of 
the elongated NH 3 distribution of L1521D (Fig. 1). This 
is consistent with the work of Goodman et al. ( 1993| ), who 
studied the occurrence of velocity gradients in a sample of 
43 ammonia cores, suggesting that the motion connected 
with the gradient, either rotation or shear, is not the ma¬ 
jor factor in the core dynamics. 

5.2. IRAS counterparts 

It is well known that the IRAS PSC offers a precious op¬ 
portunity to identify star forming regions in optically ob¬ 
scured objects, such as molecular cores. Using the IRAS 
PSC, Kenyon et al ( 1990D have recently performed a sur¬ 
vey of star forming regions in the TMC, identifying about 
100 YSOs. They have shown that YSOs are scattered in 
a large area of 15 pc x 20 pc and that star formation 
in the TMC is restricted to the dense molecular cores. 
With the aim to collect more information about the na¬ 
ture of the four molecular cores of the present list, we 
have cr oss-c orrelated their positions with the IRAS PSC. 
Clark ( 1987 ) has analysed the spatial distribution of 396 
IRAS PS located up to 2 pc from 60 molecular cores. He 
found that the distribution is strongly peaked on the am¬ 
monia cores, with the central peak reduced by an order 
of magnitude at 0.18 pc. At distances larger than this 
value, a long tail is present and the population of IRAS 
PS can be significantly affected by background sources. 
Considering that the typical ammonia core diameter is 
0.1 pc, 0.18 pc corresponds to 3.6 radii. Taking into ac¬ 
count this result, we have looked for a coincidence within 
a radius of 1.8 times the size of the four molecular cores 
of the present sample. Moreover, IRAS sources known to 
be associated with extragalactic sources have been ex¬ 
cluded (IRAS 1985| ). The result of the cross-correlation 
shows that L1521F and L1507A have no IRAS coun¬ 
terparts within the selected radius. We find two IRAS 
PS around L1521D, IRAS 04181+2655 (at a distance of 
108") and IRAS 04181+2654 (138"), and one near L1524, 
IRAS 04263+2426 (120"). It is worth noting that all these 
IRAS PS satisfy the typical IRAS colour distribution for 
IRAS counterparts of molecular cores, since they have [25- 


12] ([i-j] = log [F\/Fj], where Fy are the IRAS flux densi¬ 
ties in bands i and j, in /im) in the range between 0.42 and 
0.47 and [60-12] between 0.60 and 1.23 (for comparison, 
see figure 3 of Codella & Palla 1995). 


We thus conclude that L1521D and L1524 are associ¬ 
ated with at least one IRAS counterpart (and therefore 
presumably with YSOs), while L1521F and L1507A ap¬ 
pear not to be connected with any IRAS PS. Following 
Henning et al. ( 1990| ), the IRAS luminosities of the asso¬ 
ciated IRAS PS has been computed. The values of Tfir 
are 0.3 Lq for IRAS 04181+2655 and IRAS 04181+2654 
(L1521D) and 6.1 L 0 for IRAS 04263+2426 (L1524). This 
confirms that the ammonia molecular cores are sites of 
low-mass star formation, in accord with our results from 
the ammonia linewidths (Sect. 5.1). 


5.3. NH 3 core stability 

In order to obtain estimates of the energy terms of the 
NH 3 cores of our sample, the mean NH 3 total column 
density has been obtained as the average of all observed 
positions inside the FWHP contour of the ammonia cores 
(Figs. 1-4). This allows to calculate, using the ammonia 
abundance and assuming spherical cloud geometry, esti¬ 
mates of the mean hydrogen density and of the mass M 
(M 0 ) of a molecular core (Table 6 ). It is worth noting 
that, since M has been derived from the NH 3 total col¬ 
umn density without virial equilibrium, we have made no 
assumption about the stability of the NH 3 cores. Ammo¬ 
nia abundances in interstellar molecular clouds cover a 
wide range: the [NH 3 ]/[H 2 ] ratio is believed to vary from 
a few 10 -8 in small dark clouds up t o 10~ 5 in the dense 
cores like Orion-KL (e.g. Ho & Townes 1983| an d references 
therein). In particular, Benson & Myers ( 1983| ) studied the 
NH 3 abundance in a sample of quiescent dense molecular 
clouds and found that [NH 3 ]/[H 2 ] ranges between 3 10 -8 
and 2 10 -7 . 

In order to derive the parameters reported in this pa¬ 
per, we have assumed [NH 3 ]/[H 2 ] = 10 -7 . The resulting 
masses lie between 0.2 and 1.0 Mq , placing the NH 3 cores 
with the lowest M values (L1507A and L1521D) towards 
the edge of the mass distribution for ammonia cores in the 
TMC, which ranges from a fraction to a few tens of solar 
masses, with a median value of 4 Mq (Benson & Myers 


1989|) . 


Studying CsH 2 and H 13 CO + in TMC molecular cores, 
Mizuno et al. (1994) found that, contrary to the NH 3 re¬ 


sults of Benson & Myers ( 1989| ), the cores without young 
stars (i.e. IRAS counterparts) are less dense than those 
with stars. T he res ults of Mizuno et al. ( 1994| ) and Ben¬ 
son & Myers ( 1989 ) are consistent (average values) regard¬ 
ing the hydrogen density of the cores not connected with 
stars (4 10 4 and 3 10 4 cm -3 , resp ective ly), while for those 
which host YSOs Mizuno et al. ( 1994 ) find nn 2 — 5 10 5 


i—3 


cm °, a value larger than that given by Benson & My- 
1 10 4 cm -3 . It is worth noting that Mizuno et 


ers 


(£H). 
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Table 5. Derived parameters from NH 3 observations 


Name 

T ro t 

(K) 

T 

ex 

(K) 

JVn 

( 10 14 cm" 2 ) 

Wot 

( 10 14 cm” 2 ) 

Av 

(km s 1 ) 

(km s -1 ) 

At?NT 

(km s -1 ) 

L1521D 

9.9 (1.2) 

6.0 ( 0 . 2 ) 

1.2 (0.3) 

5.4 (1.8) 

0.27 (0.02) 

0.16 ( 0 . 02 ) 

0.22 ( 0 . 02 ) 

L1521F 

9.1 (1.0) 

9.0 (0.3) 

3.8 (0.3) 

20.1 ( 1 . 6 ) 

0.17 (0.01) 

0.16 ( 0 . 01 ) 

0.06 ( 0 . 01 ) 

L1524 

7.9 (0.7) 

7.5 (0.2) 

3.8 (0.3) 

27.7 ( 6 . 6 ) 

0.18 ( 0 . 01 ) 

0.15 (0.01) 

0.10 ( 0 . 01 ) 

L1507A 

8.9 (1.7) 

6.8 (0.3) 

1.5 (0.2) 

8.3 (3.6) 

0.20 ( 0 . 01 ) 

0.15 (0.02) 

0.13 (0.02) 


al. ( 1994| ) determined the hydrogen density from H 13 CO + 
data, while the Benson & Myers (1989) analysis is based 


on a large ve locity gradient model for the NH 3 cores. 
Mizuno et al. ( 1994 ) suggest that the disag reem ent of their 
results with those of Benson & Myers ( 1989| ) could be 
caused by differences in angular resolutions: 20 " vs. 1 '. 2 , 
respectively. We stress that our results have been obtained 
with a resoluti on of 4 0", which is still lower that that of the 
Mizuno et al. ( 1994| ) observations, but whi ch is definitely 
higher than that used by Benson & Myers ( |l989 ). The hy¬ 
drogen densities derived from our data show no remark¬ 
able difference between the 71 h 2 values of NH 3 cores with 
(L1521D and L1524) or without (L1521F and L1507A) 
IRAS counterparts. 

In Table 6 we display the energy terms of the cores of 
our sample: turbulent (.Etui J)> thermal (E t ^; J) and grav¬ 
itational (E gr ; J). The last column is for the parameter a 
= E gr /(E tu + Eth), which is a measure of core stability. 
The turbulent, thermal and gravitational (assuming a ho¬ 
mogeneous sphere) energy terms have been derived using 
the formulae reported by Harju et al. ( 1993 ). 

Star formation occurs most readily in molecular 
clumps in which the internal motions are sufficiently re¬ 
duced to allow the occurrence of gravitational collapse. 
Table 6 shows that the thermal energy is always larger 
than the turbulent one. The values of a for the four cores 
lie between 0.2 (L1521D) and 1.2 (L1524); a = 1 implies 
gravitational equilibrium in the absence of a magnetic 
field. The values obtained here are only slightly differ¬ 
ent from unity. Taking into account the contribution of 
a magnetic field to the energetic support against gravi¬ 
tational collapse, considering also uncertainties in the as¬ 
sumed [NH 3 ]/[H 2 ] value and deviations from the adopted 
(spherical) geometry of the cores, we can conclude that the 
four ammonia cores are close to equilibrium. These con¬ 
siderations are i n agreement with t he re sults of Myers & 
Goodman ( 198S ) and Harju et al. ( 1993 ) who found that 
the NH 3 cores in the TMC are approximately virialized. 


5-4■ Derived parameters from HC 5 N observations 

An estimate of the total HC 5 N column density (-/Vhc 5 n; 
cm -2 ) has been derived through the standard equation 


(e.g. Olano et al. |1988|) , using t he m olecular parameters 
give n, e.g ., by Alexander et al. ( 1976 ) and Churchwell et 
al. ( 1978|) . Table 4 lists the derived cyanodiacetylene col¬ 
umn densities at the ammonia peak positions of the four 


molecular cores: the values range from 1.6 10 12 to 9.2 10 12 
cm~ 2 , in agreement with the -/Vhc 5 n values of othe r mol ec- 
ular cores located in the TMC (Benson & Myers 1983 ). 

In order to estimate the HC 5 N abundances a compar¬ 
ison between HC 5 N and NH 3 column densities has been 
made. We presently do not know the size and geometry 
of the HC 5 N distribution in the ammonia cores. However, 
the HC 5 N and NH 3 maps given by Benson & Myers (1983) 
show that the size and shape of the maps are similar and 
that the peak positions differ significantly only for the 
TMC-1 molecular core. Therefore, we have estimated the 
[HCsNJ/pd^] ratio using the cyanodiacetylene and ammo¬ 
nia column densities and the NH 3 abundance. The ob¬ 
tained [HC. 5 N]/[NH 3 ] values are: 8.3 10 -3 (L1521D), 4.6 
1CT 3 (L1521F), 2.7 1CT 3 (L1524) and 10.7 10 ~ 3 (L1507A). 
Thus, the derived HC 5 N abundances relative to H 2 are: 8.3 
1CT 10 (L1521D), 4.6 10 - 10 (L1521F), 2.7 10 " 10 (L1524) 
and 1.1 10 -9 (L1507 A), in agreement with the suggestion 
of Benson & Myers ( 1983 ), that relative cyanodiacetylene 
abundances in dense cool cores vary from ~ 10~ 10 to ~ 
1 CT 9 . 


6. Conclusions 

We have discovered 4 dense cores (L1521D, L1521F, L1524 
and L1507A) in the TMC, increasing the number of known 
such objects to 20 in this prototypical dark cloud complex. 
The main results of our study can be summarized as fol¬ 
lows: 

1. The ammonia observations reveal clumps with sizes 
of 0.06-0.09 pc and yield kinetic temperatures be¬ 
tween 7.9 and 9.9 K. The actual intrinsic linewidths 
are between 0.17 and 0.27 km s , while the thermal 
linewidths are around 0.15 km s _1 . Therefore, the non- 
thermal velocity component is rather small. The ratio 
between the intrinsic and thermal linewidth is about 
1.3, indicating that only low-mass stars can eventually 
form in these NH 3 cores. This is also supported by the 
association of two of the NH 3 cores with low luminosity 
IRAS counterparts (Tfir < 10 Eq). 

2. The NH 3 data show that the measured (1,1) excitation 
temperatures are only slightly smaller than the rota¬ 
tional temperatures suggesting that the ( 1 , 1 ) inversion 
transition is close to being thermalised. From the peak 
NH 3 emission hydrogen densities (between 0.6 and 19.9 
10 4 cm -3 ) have been estimated. Moreover, from the 
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Table 6. Derived parameters related to the core stability of the NH3 condensations 


Name 

n H 2 

(10 4 cm" 3 ) 

M 

(M Q ) 

Etu 

(10 34 J) 

Eth 

(10 34 J) 

Egr 

(10 34 J) 

a 

L1521D 

0.6 

0.3 

0.9 

3.5 

0.9 

0.2 

L1521F 

19.9 

0.7 

0.8 

7.9 

7.4 

0.9 

LI 524 

2.5 

1.0 

1.0 

8.9 

11.8 

1.2 

L1507A 

1.1 

0.2 

0.2 

2.2 

1.2 

0.5 


average of all observed positions inside the FWHP con¬ 
tour of the ammonia cores, the core masses ( 0.2 - 1.0 
Mq) have been estimated. The derived hydrogen den¬ 
sities show no remarkable difference between NH 3 cores 
with or without IRAS counterparts. The low M values 
place the four objects towards the edge of the M dis¬ 
tribution for ammonia cores in the TMC, which ranges 
from fractions to few tens of solar masses, with a me¬ 
dian value of 4 Mq. 

3. Estimates of the turbulent, thermal and gravitational 
energy have been derived. The thermal energy is al¬ 
ways larger than the turbulent one. The four cores may 
be close to equilibrium. 

4. Emission at 24.0 GHz, caused by the J=9-8 transition 
of HC 5 N, has been detected towards the ammonia peak 
positions of the four molecular cores. This increases the 
number of known such objects from 11 to 15 in the 
TMC (cf. Benson & Myers 1980, 1983). The derived 
HC 5 N column densities range between 1.6 10 12 and 
9.2 10 12 cm" 2 , in agreement with results of previous 
HC 5 N observations in the TMC. 
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